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Abstract 

We investigate the idea of a "general boundary" formulation of quantum field theory in the context of 
j^. , the Euclidean free scalar field. We propose a precise definition for an evolution kernel that propagates 

£NJ ■ the field through arbitrary spacetime regions. We show that this kernel satisfies an evolution equation 

which governs its dependence on deformations of the boundary surface and generalizes the ordinary 
(Euclidean) Schrodinger equation. We also derive the classical counterpart of this equation, which is a 
, Hamilton- Jacobi equation for general boundary surfaces. 
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rn ' 1 Introduction 

O ■ 

In quantum field theory (QFT) on Minkowski space, we can use the Schrodinger picture and have states 
"V* ' associated to flat spacelike (hyper-)surfaces. The transition amplitude between an initial state and a final 
^■i' state is obtained by acting with the unitary evolution operator on the former and taking the inner product 
with the latter. The possibility of a Schrodinger picture has also been considered in QFT on curved spacetime 
[1]. In this case, states live on arbitrary spacelike Cauchy surfaces forming a foliation of spacetime. Evolution 
along these surfaces is non-unitary in general, as it does not correspond to a symmetry of the metric. In 
background independent quantum gravity, on the other hand, there is no fixed spacetime geometry; in this 
case, states live on arbitrary Cauchy surfaces and the requirement that the surface is spacelike is encoded 
in the state itself, which represents a quantum state of a spacelike geometry (see e.g. [2, 3]). 

In all these cases, transition amplitudes are calculated for boundary states (i.e. an initial and final 
state) defined on spacelike boundaries. Recently, Oeckl has suggested that it could be possible to relax 
this restriction to spacelike boundaries in QFT [4, 5]. Oeckl offers heuristic arguments which suggest that 
transition amplitudes can be associated to a wider class of boundaries, as we do in topological quantum 
field theory [6]. These more general boundaries may include hypersurfaces which are partially timelike, 
that enclose a finite region of spacetime, or disjoint unions of such sets. This would imply, for instance, 
that in theories like QED or QCD, we could associate quantum "states" to a hyper sphere, a hypercube or 
more exotic surfaces, and assign probability amplitudes to them. Similar suggestions were made in [3], with 
different motivations. 

This "general boundary" approach to QFT could be interesting for several reasons. Firstly, finite 
closed boundaries represent the way real experiments are set up more directly than constant-time surfaces. 
A realistic experiment is confined to a finite region of spacetime. In particle colliders, for instance, the 
interaction region is enclosed by a finite outer region where state preparation and measurement take place. 
As sketched in Fig. 1, the walls and openings of a particle detector trace out a hypercube in spacetime. A 
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Figure 1: Spacetime diagram of particle scattering. 

"state" on the hypercube's surface would represent both incoming beams and jets of outgoing particles in 
a completely local fashion, without making any reference to inaccessible infinitely distant regions. 

Secondly, in a quantum theory of gravity closed boundaries may provide a way to define scattering 
amplitudes, and help in solving the traditional interpretational difficulties of background independent QFT. 
This idea has been recently studied in [7], where it has also been used to propose an explicit way for 
computing the Minkowski vacuum state from a spinfoam model [8]-[ll]. In a background independent 
theory the conventional spacelike states do not impose any constraint on the proper time lapsed between 
the initial and final states. As a result, the transition amplitude stems from a superposition of processes 
whose duration may range from microscopic to cosmic time scales. Fixing a timelike boundary can control 
the time lapsed during the experiment. 

Furthermore, in spinfoam approaches, the introduction of general boundaries might open up the pos- 
sibility of quantizing 3-geometries along time-like surfaces and clarify the physical meaning of Lorentzian 
spinfoams. 

Finally, a general boundary formulation could give us a broader perspective on QFT: it would stress 
geometrical aspects of QFT by no longer singling out a special subclass of surfaces, and may shed some 
light on the holographic principle, which states that the complete information about a spacetime region can 
be encoded in its boundary. 

As noted by Oeckl, a heuristic idea for adapting QFT to general boundaries is provided by Feynman's 
sum-over-paths-picture. Given an arbitrary spacetime region V, bounded by a 3d hypersurface, the Feynman 
path integral over the spacetime region V, with fixed boundary value (p of the field, defines a functional 
W[ip, V]. This functional can be seen as a generalized evolution kernel, or a generalized field propagator. 
The path integral is therefore a natural starting point for developing a general boundary formalism. 

The path to make these ideas precise is long. There are two types of problems. Firstly, the probabilistic 
interpretation of quantum theory and QFT must be adapted to this more general case. The physical meaning 
of states at fixed time and their relation to physical measurements are well established; the extension to 
arbitrary boundaries is probably doable, but far from obvious. It requires us to treat quantum state 
preparation and quantum measurement on the same ground, and to give a precise interpretation to the 
general probability amplitudes. Some steps in this direction can be found in [5] and [3]. 

Secondly, the mathematical apparatus of QFT, i.e. the path integral and operator formalism, needs to 
be extended to general spacetime regions. On a formal level, such a generalization appears natural for path 
integrals, but it is far from clear that it can be given a concrete and well-defined meaning. 

In this paper we focus on the second of these issues: the definition of the field theoretical functional 
integral over an arbitrary region, and its relation to operator equations. We start to address the problem 
by considering the simplest system: Euclidean free scalar field theory. In this context, we propose an 
exact definition for the propagator kernel W[ip, V], based on limits of lattice path integrals. Under a 
number of assumptions, we can show that the propagator satisfies a generalized Schrodinger equation, 
of the Tomonaga-Schwinger kind [12, 13]. The equation governs the way the propagator changes under 
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infinitesimal deformations of V. It reduces to the ordinary Schrodinger equation in the case in which a 
boundary surface of V is a constant-time surface and the deformation is a global shift in time. 

With this result, we provide a first step towards constructing an operator formalism for general bound- 
aries. The derivation can be seen as a higher-dimensional generalization of Feynman's path integral deriva- 
tion of the Schrodinger equation for a single particle [14]. The main assumption we need is the existence of 
a rotationally invariant continuum limit. 

We also derive the classical counterpart of the evolution equation: a generalized version of the Euclidean 
Hamilton-Jacobi equation. At present, we have no prescription for Wick rotation, so we cannot give any 
Lorentzian form for the propagator or the evolution equation. Hints in this direction were given in [7]. 

If one continues along this line, the ultimate goal would be to construct a full general boundary formal- 
ism for background dependent QFTs, which incorporates Wick rotation, interactions and renormalization. 
While of interest in itself, such a project could be also viewed as a testing ground for the general boundary 
method, which would prepare us for applying it in the more difficult context of background free QFT: there, 
as indicated before, the use of generalized boundary conditions may not only be helpful, but also essential 
for understanding the theory. 

Our technique for deriving the evolution equation could be of interest in view of the attempts to relate 
canonical and path integral formulations of quantum gravity, i.e. when deriving the Wheeler-DeWitt equa- 
tion from a concrete realization of a sum over geometries. (For existing results on this problem, see e.g. 
[15].) 

The paper is organized as follows. In section 2, we present some of the heuristic considerations about 
state functionals on general boundaries, and their associated evolution kernel. Section 3 deals with the 
classical case: we present two derivations of the generalized Hamilton-Jacobi equation. The lattice regu- 
larization of the quantum propagator is defined in section 4. There, we also state the assumptions which 
are then used in section 5 for deriving the generalized Schrodinger equation. Both Hamilton-Jacobi and 
Schodinger equation are given in their integral form. In the appendix we clarify the relation with the local 
notation in [3]. 

Notation. V is the spacetime domain over which the action and the path integrals are defined. E is 
the boundary of V. The letter cf> denotes a real scalar field on V, while ip stands for its restriction to E, 
i.e. <p = (f)\-£. Depending on the context, <p can be a solution of the classical equations of motion or an 
arbitrary field configuration. The action associated to <j> is written as S[(f>, V]. When boundary conditions 
(cp, E) determine a classical solution cf> on V, we denote the corresponding value of the action by S[tp, V\. 
Thus, the functional S[ip, V] can be viewed as a Hamilton function (see sec. 3.3 of [3]). Vector components 
carry greek indices fi,v,... (e.g. v = (v^)). The dimension of spacetime is d. The symbol J v d d x represents 
integrals over V, while integrals over E are indicated by J s dE(x). The letter n denotes the outward pointing 
and unit normal vector of E. The normal derivative is written as d n , while Vs is the gradient along E. 
Accordingly, the full gradient V decomposes on E as 

V| s = n<9 ri + V s . (1) 

In section 4, we introduce a lattice with lattice spacing a and regularize various continuum quantities. Their 
discrete analog ues are designated by the index a: for example, ip, V and S become (p a , V a and S a . 

2 General Boundary Approach 

What is the meaning of a state on a general surface which is not necessarilly spacelike? What does it mean 
to propagate fields along a general spacetime domain? Following [5], an intuitive answer to these questions 
is provided by the path integral approach to QFT. We illustrate this intuitive idea in this section, as a 
heuristic motivation for the more rigorous definitions and developments in the remainder of the article. For 



3 



Figure 2: From V/j to general V. 



simplicity, we refer here to a scalar field theory, but similar considerations can be extended to various path 
integral formulations of QFT, including sum over metrics or spinfoam models in quantum gravity. 

Consider Minkowskian scalar QFT in the Schrodinger picture. Let be an initial state at time U 
and \^f) a final state at time tf. The transition amplitude between the two is 

A = (^ f \e- iH ^f- ti ^ h \^ i ). (2) 
Using the functional representation, the amplitude (2) can be expressed as a convolution 

a = jD<p f J zm*^/W^/'*/;^]^N ( 3 ) 

with the propagator kernel 

W^tf^U] := Me-^f-Wfa) . (4) 

This field propagator is a functional of the field: it should not be confused with the Feynman propagator, 
which is a two-point function, and propagates particles. When rewritten as a path integral, this kernel takes 
the form 

W[<p f ,tf,<p i ,t i ]= J D<j>e iS ^y n . (5) 

4>{-,ti) = tpi , 
4>i,tf) = <p f 

The action integral extends over the spacetime region Vji := R d_1 x and the path integral sums over 

all field configurations 4> on V/j that coincide with the fields ipj and (pi on the boundary. The complete 
boundary consists of two parts: the hyperplane Sj at the initial time ti, and the hyperplane at the final 
time tf. We call their union Sjj := U Sj. If we view <pf and (pi as components of a single boundary field 
if fi := (ipf,<fi) on Tifi, we can write the evolution kernel (5) more concisely as 

W[<p f i,V fi ] := / D<f> e iS ^ v f^ h . 

H\j: Si =ffi 

With this notation, it seems natural to introduce a propagator functional for more general spacetime regions 
V (see Fig. 2): we define it as 

W[p,V\ := / D<t> e iS[<t> ' V]/h . (6) 

J<l>\-Z=ip 

Here <p varies freely on the interior of V and is fixed to the value ip on the boundary S. Of course, this is 
only a formal expression, and it is not clear that it can be given mathematical meaning. Let us suppose for 
the moment that it has meaning and see what would follow from it. 

Ordinary propagators satisfy convolution (or Markov) identities which result from the subdivision or 
joining of time intervals. If the functional W behaves the way our naive picture tells us, the splitting and 
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Figure 3: Splitting of V. 
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Figure 4: Splitting of V/j. 



joining of volumes should translate into analogous convolution relations. For instance, if V is divided as 
shown Fig. 3, the new regions V cb and V ba carry propagators 

W[& c ,w),V cb ] = J D<t>e lS ^ v ^ h , (7) 

0|e c6 =(¥>o¥>&) 

W[(<p b ,<p a ),V ba ] = J D4>e iS ^ n . (8) 

4>\s ba =(<Pb,>Pa) 

When integrating the product of (7) and (8) over the field (p b along the common boundary, one recovers 
the original propagator: 

W[(<p c , <p a ), V] = J D Vb W[(<p c , Vb ), V cb ] W[(<p b , <p a ), V ba ] . (9) 

Similarly, the infinite strip V/j between tf and U could be cut by a "middle" surface S m as in Fig. 4, giving 
the new volumes Vf m and V m i- Its kernel decomposes as 

W[{<Pf, <Pi), V fi ] = J D(f m W[(<pf, <p m ), V fm ] W[(<f m , ipi), V mi ] . 

Thus, the evolution of to the final time is divided into two steps: using the propagator on V m i, we 
evolve up to the surface S m and obtain the intermediate state 

^mYPm] ■= / d(piW[(<p m ,ipi),V mi \^i[ipi\ . 
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Figure 5: Evolution to a closed surface £. 



The kernel W[.,V f m ] covers the remaining evolution and gives the original amplitude (2) when convoluted 
with $>* f and * m 

A = J D(p f J D(p m V* f [<pf]W[((pf,(p m ),Vf m ]V m [<Pm]- (10) 

Since the same amplitude can either be calculated from and or from tyf and \I/ m , the wave functional 
^> m encodes all physical information about the intial state. On account of this property, we say that ^ m 
is the state functional which results from evolving by the volume V m i. As for ordinary state functionals, 
one can think of ^ m as being an element \^ m ) in a Hilbert space, which we call Ws m - The latter consists 
of functionals of fields over S m and has the inner product 

0N*l) := / D<p m *5[Vm]*l[Vm] , |*1> , 1*2) G H Sm • 

It is important to note that the evolution map from TC^ to Ws m need not be unitary. The results of 
Torre and Varadarajan show, in fact, that in flat spacetime state evolution between curved Cauchy surfaces 
cannot be implemented unitarily [16]. Nevertheless, a probability interpretation is viable for states in H^ m , 
as the meaning of amplitudes such as (10) can be traced back to that of the standard amplitude (3). 

Consider now a more unconventional example. Cut out a bounded and simply connected set V m from 
Vfi and denote the remaining volume by V/ m j (Fig. 5). This time we define the state functional ^> m by 
"evolving" both * j and ^> m to the middle boundary S m , i.e. 

^mVfim] ■= J D(f f J Difii ^*f[(pf]W[((pf,(p m ,(pi),Vf m i\^i[(pi\. 

Clearly, the amplitude (3) is now equal to 

A = J Dif rn W[ip m , V m ] V m [<Pm} ■ (11) 

Therefore, the functional fy m contains the entire information needed to compute the transition amplitude 
between and 

To make this more concrete and more intuitive, suppose that the scalar field theory is free and that 
Vl/j and tyf are the initial and final one-particle states of a single, localized particle whose (smeared out) 
worldline passes through V m . In both functionals, the presence of the particle appears as a local deviation 
from the vacuum, in the functional dependence. Likewise, it is natural to presume that the functional form 
of * m reflects where the worldline of the particle enters and exits the volume V m . 

How can we interpret the "state" ^ m and the associated amplitude (11)? To answer this, let us get 
back to equation (2). Notice that the amplitude A depends on the couple of states |^/))- This 

couple represents a possible outcome of a measurement at time tf as well as a state preparation at time 
t, L . A state preparation is itself a quantum measurement, therefore we can say that this couple represents 
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a possible outcome of an ensemble of quantum measurements performed at times U and tf. We may 
introduce a name to denote such a couple. We call it a process, since the two states 1^/)), taken 

together, represent the ensemble of data (initial and final) that we can gather about a physical process. A 
probability amplitude is associated to the entire process (|*j), |^/))- Now, it is clear that the functional \I> m 
represents a generalization of this idea of a process. It is tempting to presume that ^> m can be interpreted 
as representing a possible outcome of quantum measurements that can be made on E m . In the example of 
the particle above, for instance, it will represent the detection of the incoming and outgoing particle. 

The idea is that given an arbitrary closed surface, the possible results of the ensemble of measurements 
that we can make on it determines a space of generalized "states" which can be associated to the surface. 
Each such state represents a process whose probabilistic amplitude is provided by expression (11). The 
conventional formalism is recovered when the surface is formed by two parallell spacelike planes. For more 
details on the physical interpretation of general boundary states, see sec. 5.3 of [3]. 



2.1 Operator Formalism 

If path integrals can be defined for general boundaries, how would a corresponding operator formalism 
look like? In particular, is there an operator that governs the dynamics, as the Hamiltonian does for rigid 
time translations? Recall that the Hamiltonian can be recovered from the path integral by considering an 
infinitesimal shift of the final time. For example, if we displace by a time interval At the final surface E f 
in (5), keeping the same boundary field iff, the new propagator results from the convolution 

W[<p f ,t f + At;tpi,ti] = J DtpW[ip f ,t f + At;tp,t f ]W[tp,t f ;ipi,ti]. (12) 

For infinitesimal At, this gives the Schrodinger equation, which expresses the variation of W in terms of 
the Hamiltonian operator 

lih^--H[ip f ,-ih— ]) W[<p f ,tf,<pi,ti]=0, (13) 



dt f ; " diff 



where 



H[<p f ,* f ] = £ dS 1 (ttJ + (V m ) 2 + mV/) • (14) 

Similarly, if iff is displaced in a tangential direction ey along Ej, the variation of W is generated by the 
momentum operator 

e|| . P[<^ , 

where 



P b/> 7T/] = - / dY,(x) V^f(x) TT f . 



In the case of a general volume V, it is natural to expect that deformations of the boundary surface E lead 
to an analogous functional differential equation for the propagator. However, for a general shape of V there 
is no notion of preferred rigid displacement of the boundary. We must consider arbitrary deformations of the 
boundary surface, and we expect that the associated change in W is governed by a generalized Schrodinger 
equation (see Fig. 2). In the same way that H and P generate temporal and spatial shifts, the operators in 
such a Schrodinger equation could be seen as the generators for general boundary deformations of W. In a 
difformopshim invariant QFT, the analogous W- functional would be independent of E, and the generalized 
Schrodinger equation reduces to the Wheeler-DeWitt equation. 

In the present paper, we consider only Euclidean field theory, so we seek to define the Euclidean form 

W[<p,V] := / D<f> e- s ^' n . (15) 
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of the propagator (6), and generalize the Euclidean version of the Schrodinger equation (13). 

Before dealing with path integrals and deformations of their boundaries, however, we discuss the 
analagous problem in classical field theory. The classical counterpart of the Schrodinger equation is the 
Hamilton-Jacobi equation. The Hamilton function S[ipf, tf, ipi, ij] is a function of the same arguments as 
the field propagator (4). It is defined as the value of the action of the classical field configuration which 
solves the equations of motion and has the given boundary values. It satisfies the Hamilton-Jacobi equation 

d & S 

—S[cp f ,t f ;i Pi ,t i } + H[ip f ,—}=0. (16) 

For more general regions V, the Hamilton function becomes a functional of V and the boundary field ip 
specified on E, In the next section we show that this functional satisfies a generalized Hamilton-Jacobi 
equation which governs its dependence on arbitrary variations of V. 



3 Generalized Hamilton-Jacobi Equation 

Let V be an open and simply connected subset of Euclidean d-dimensional space R, d . We consider the 
Euclidean action 

- ! 1 

(17) 



S[<f>, V]= f d d x 
Jv 



I(V0) 2 + ImV + f/(^) 



where U is some polynomial potential in (f>. In the classical case, unlike in the quantum case, an interaction 
term can be added without complicating the derivation that follows. The equations of motion are 

BU 

n<f ) -m 2 ( f>- — = 0. (18) 

The Hamilton function S[(p, V] is defined by S[ip, V] = S[(f),V], where 4> solves (18) and (f>\% = ip. It is 
defined for all values (ip, S) where this solution exists and is multivalued if this solution is not unique. 

We now study the change in S[<p, V] under a local variation of V. To make this precise, consider a 
vector field N = (iV 7 *) over R d . N induces a flow on R d which we denote by a : R x R d — ► 1R, d . Define the 
transformed volume as V s := a(s, V). Likewise, S s := a(s, £). 

To define the change in S[np, V] under a variation of V, we need to specify what value the boundary 
field should take on the new boundary S s . We choose it to be the pull-forward by a s = a(s,.), i.e. 
ip s := ip o cj^ 1 = (T s ^ip. Let us assume that the point (ip, S) is regular in the space of boundary conditions, in 
the sense that slightly deformed boundary conditions (99 s , X s ) give a new unique solution s on V s , close to 
the previous one. In this case, the number S[ip s , V s ] is well-defined and we can write down the differential 
quotient 

L N S[ip, V] := lim - (S[ip s , V s ] - S[<p, V]) , (19) 

s^O S 

with the vector field N as a parameter. As we show below, this limit exists and the map is a functional 
differential operator. The local form of this differential operator is given in the appendix. 
We decompose the restriction of N to E into its components normal and tangential to S, 

7V |S = N±n + TV,, , 

where the scalar field iV_|_ is defined as 

N ± := . 

Observe that under a small variation 5ip of the boundary field, we have 

SS[tp, V] = S[Scp, V] = S[5<t>, V] 



-L 



(fx 

v 



dU 

V^Scf) Vu<p + m 2 (p 5cj) + —5(f) 

Oq) 
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V s 



Figure 6: Definition of p^. 



= [ d d 
Jv 



Therefore, we have 



3.1 Direct Derivation 



cffi d n 4> 54> = 
5S 



5<p{x) 



dTi d n 4> Sp . 
[p, V] = d n (f>(x) . 



(20) 



Suppose for the moment that V is only extended by the deformation (i.e. V C V s for every s). Then, the 
most direct derivation of the Hamilton-Jacobi equation can be obtained by considering the restriction p^ 
of (f) s to S: that is, the value of the classical solution on £ when the boundary condition p s is specified on 
S s (see Fig. 6). 

Note that ip^ = <p° = <p. By inserting S[ip^, V] into the difference, i.e. 

S[p s , V s ] - S[p, V] = S[p s , V s ] - S[&, V] + S[p s s , V] - S[p, V] , 
the differential quotient becomes a sum of two limits: 

Hm S[p s ,V s ]-S[p,V] = ^ S[<p*,V] - S^V] + ^ S[p° s ,V}-S[p,V} 

s— >o S s^O s s^O s 

As p s and <p^ are part of the same solution, the first limit is easily seen to be 



lim 



S[^,V s ]-S[p^,V] 



In the last line, we used the decomposition (1) and equation (20). The second differential quotient gives 

S[p^,V]-S[p,V] 



lim 



f 6S d s 1 
= J^dZ — (-N ± d n <f> - N\\ ■ Vxp) . 



Altogether one has 

L N S[<p,V]= J^dE^ 



N ± 



(21) 
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We started from the assumption that V C V s for all s, but it is easy to see that the previous argument can 
be adapted to the general case where the volume V is partly extended and partly decreased. 
If we introduce the quantities 

H N [<p, tt, V] := ^ dS Nj_ (-^vr 2 + ^(V s ^) 2 + ^mV + U(<p)J , 

P N [ip,TT,V] := -Jd£N\\ -Vs^tt, 

equation (21) takes the form 

L N S[p, V] = H N [<p, S ^V]+P N [ V ,^,V}. (22) 

When V is a strip of spacetime between times U and tf, and N±\t f = 1, N±\ ti = 0, iVy = 0, equation (22) 
reduces to the usual Hamilton-Jacobi equation 

■^S[<pf,t f ;<p i ,t i ]= J^dZ +\{^}? + \m\) + U^ f ) 

of Euclidean field theory. Hence (22) can be seen as a geometric generalization of the Hamilton-Jacobi 
equation. 

3.2 Alternative Derivation 

Let us describe another way of evaluating the "deformation derivative" L]yS[(p,V]. The spacetime metric 
tensor g enters in the definition of the action and therefore in the definition of S[(p, V\. Let us write this 
dependence explicitely as S[(p,g,V]. A diffeomorphism that acts on 4>, the boundary S and the metric g, 
leaves the action invariant. Therefore 

S[ip s ,g s ,V s ]=S[ip,g,V]. 

Equivalently, 

S[<p>,g,V a ]=S[<p,g-,V]. 
Plugging this into the definition of the operator (19) gives 

L N S[<p, g, V] = km - (S[<p, <T S , V] - S[<p, g, V}) , (23) 

s— >0 S 

which is a variation of the action w.r.t. the metric only. Now we can use the definition of the energy- 
momentum tensor to obtain 

L N S[cp,g,V] = I J v d d xT^^g-°\ s= o = ^ J v d d xT^(-V^N u -V u N^) 

= -J d d xT^ u V^N u = - J dEn li T' u 'N v + J d d xV ll T^N v 

=o 

= -J^d d xn^ v N u (24) 

In the last two steps we used Stoke's theorem and the equations of motion respectively. On the other hand, 
we know that 
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and 



-N ± 
-N ± 



-\(d n 4>? + ^(V s 0) 2 + + U{4>) 



+ d n 0(N ± d n <f> + TV,, • V s 0) 



Inserting (25) in (24) and using (20), we arrive again at the generalized Hamilton- Jacobi equation 



L N S\f, V] = / dS I N ± 



■iQ ,+ W> ,+ ^ +i "« 



" dip 



(25) 



4 Definition of the Evolution Kernel 

In this section, we define a Euclidean free field propagator for arbitrary spacetime domains V. Limits of 
lattice path integrals are used to give a precise meaning to the expression (15). We begin by considering 
the case V = Vfi and derive the lattice path integral from the operator formalism. Then, we propose a way 
to extend this expression to more general volumes V. 



4.1 Prom Operators to Path Integrals 

The transition from operator formalism to path integral is a standard procedure. We repeat it here, since 
treatments of lattice field theory usually omit normalization factors. There, constant factors drop out when 
dividing by the partition function Z. In our case, their precise form will be crucial for the definition of the 
propagator. 

In the Schrodinger picture, the space of states H is associated to the manifold R d-1 : we regularize it 
by a finite lattice 

S a := {x G aZ d_1 | - Ma < \xi\ < Ma , i = 1, ... ,d - 1} 

with lattice constant a > and edge length 2aM, MeM. ej is the unit vector in the ith direction. For a 
scalar function / on S a , the forward derivative is 



Vi/(f) := 



4>(x + aej) - 4>(x) 



and we set 4>{x + aej) := <f>(x — aMei) when Xi = aM. Let {4>(x)}, {n(x)} be canonical operators with 
eigenstates {| </>)}, {K)} such that 



and 



4>(x)\<t>) =<f>(x)\<f>), tt(x)\tt) = ir(x)\ir) , x €. S a . 



ih 



[4>{x)^{y)] = -jztS(x- y) , x,y£S a . 
a a 1 

The eigenstates are normalized as 

{W) = n Km - cp'(x)) , (7r,7r') = I] - vr'(x)) , 

and give rise to completeness relations 



(26) 
(27) 

(28) 



n 

y ses a J -' x 



d<f>(x) \ 1^(^1 = 1, 



n 



dir(x) \tt)(tt\ = 1 • 



(29) 
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From (26), (27) and (28), it follows that 



and 



The Hamiltonian operator is 



X&Sa 



U\S) + l -(y^{x) + \m^\x) 



We rewrite the Euclidean propagator 

(iffle-^f-^lifi), t f -U = na, 
by inserting repeatedly the completeness relations (29): 

we-™*^) = ( nn / ^(^l fffn / >< 

\fc=l x J / Vfe=0 x / 

x (^ / |7r„_ 1 )(7r n _ 1 |e- a// / ?l |0 n _ 1 )(0 n _ 1 |7r„_ 2 )(7r n _ 2 |e- aH / ?i |0 n _ 2 ) • • • {0! |7T ) <tt | e" a ^/ a | ^) 

After making the replacement 

e -aH/h = e ~aT/h e -aV/h + 0( Q 2) _^ g-aT/^g-aV/ft 

and using (31), we obtain 

M| II / Whiz) Ml II / dlT k{x) ) (<Pf\n n -l){Xn-l\<l>n-l)(<f>n-l\nn-2)(nn-2\<f>n-2) ■ ■ ■ 
\fc=l x / \fc=0 x ^ / 



n-l 



■ ■ ■ (0iko)(vr o |^) exp I - aH[4> k ,TT k ] 



k=0 



<t>0 = ¥>z 



nn / d^w) Pnn / (n 

U=l x J / \fe=0 x / \ x 



2ttH 



2n 



n—1 



ex Pi^E a E a(i 1 



k=0 x 



.4>k+i{x) - (t>k{x) 



Kk{x) - ^ 2 k (x) - ]^{V<p k ) 2 {x) - -m 2 4> 2 k (x) 



4>o 



Integration over the momenta yields the path integral 



\k=l x 



n-l 



X 



k=0 



1 (4>k+l{x) - (p k (x)\ 2 1 2 1 2 2 



x i 



<t>0 = <Pi , 
4>n = <fif 
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Figure 7: Lattice diagram for path integral on V/j. 



In the zeroth and nth layer, <p is fixed to the initial and final values ipi and iff respectively, while it is 
integrated over from layer 1 to n — 1, weighted by the exponentiated action. 

We can make this formula more symmetric with respect to the boundaries U and tf, if we add potential 
terms to the nth layer, writing 



o \ n/2N 



\k=l x 



,d-2 



2irh 



x exp 



n-l 

£5 



fe=0 



1 f 4>k+i(x) - <p k {x) 

2 V a 



-W 2 



fc=0 



EUC^*) 2 (*) + «m 2 *!(a 



00 = ¥>z i 
0n = V/ 



• (33) 



Clearly, such a change does not affect the continuum limit. We also rewrite the normalization factors: in 
(33), there are n factors of 



~.d— 2 



2vr?i 



(34) 



for each x £ S a . We express this in a more geometric fashion by attributing a factor C a to every spacetime 
point x = (x,ti + ka) for which 4>k{x) is integrated over, and by associating a factor \[C~ a to each point in 
the initial and final layer: 



w a [<pf,tf,<pi,u] ■■= f n f n n / ^^(5)] x 



x exp ■ 



a 



x£S a lk=0 



+ E(l^*) 2 ( £ ) + 5 w » 2 ^ 



fc=0 



00 = > 
0n = ff 



.(35) 



In Fig. 7, this is represented diagrammatically for the case d = 2: open points stand for an integration 
over the associated field variable and carry a factor C a . Boundary points are solid and contribute a factor 
\[C~ a . For each point there is a mass term in the action, and each link between points gives a term with the 
corresponding lattice derivative. The dual lattice is drawn shaded. 



4.2 General Definition 

By applying the same rules to more complicated arrangements of points, we can define a path integral 
regularization for general volumes V. Let V C R d be open and its boundary E piecewise smooth. We use 
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Figure 8: Lattice diagram for a general volume V. 



hyper cubic lattices 

L a := {x e a% d | - aM < \x^\ < aM , /j, = 1, . . . ,d} 

with lattice constant a > and edge length 2aM, M £ N. is the unit vector in the fith direction. A 
lattice point x and a direction fx define a link 

/ = (x,n) . 

The associated lattice gradient is 

V</ = V,/(z) :=fa^. 

a 

Given a subset P C L a , l(P) denotes the set of links that connect points within P. Let 

V a :=L a nV 

be the intersection of V with the lattice. The points of V a fall into three categories (Fig. 8): we call a point 
interior if it has 2d links to points of V a . If a point is linked to less than 2d points of V a , but connected to 
at least one interior point, it is a boundary point. The remaining points of V a have only links to boundary 
points and we will not use them when representing the path integral on the lattice (they are drawn shaded 
in Fig. 8). The set of relevant points is therefore 

V a := la U S a , 

where I a and S a denote the set of interior and boundary points respectively. 

On the lattice, the path integral becomes a summation over scalar fields <f> : V a — > R on V a . The action 
consists of contributions from links in l(V a ) and points in V a : 

S[cP,V a ]:= J2 «^Wf+E a 1 m ^) 

lel(V a ) X&Va 

Given a continuous boundary field (p on S, one has to translate it into boundary data for V a . We do so by 
defining the discrete boundary field 

(p a : S a -> R, ^ a (x) = (p(pmdx(x)) , 

The function pmd^ (pmd stands for "point of minimal distance") returns a point on S which has minimal 
distance to x. Now, we have all the necessary notation to give the regularized form of the propagator 
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W[ip, V]: we specify it as 

W a [p a ,V a ] := I J] VC~a) [ / II C a d4{x)\ exp (~S[tf>,V a ]) , (36) 

\xes„ / V ^e/ a J \ n j 4>\v a =<p a 

with factors C a as in (34). The continuum propagator W[(p, V] is then defined by the limit of vanishing 
lattice constant and infinite lattice size: 

W[p,V\ := lim lim W a [<p a ,V a \. 
To simplify notation, we omit the lim -symbol in the remainder of the text. That is, the limit of infinite 

M^oo 

lattice size (for constant a) is implicit in all subsequent formulas. 

We now make a number of unproven assumptions about the regularization (36): 

(Al) The propagator (36) has a continuum limit: that is, there is a non-trivial space -F(S) of boundary 
fields on £ such that for each ip G the limit 



W[<p,V] := limW a [<p a ,V a ] 

a— »0 



is well-defined. 



(A2) W reproduces the conventional propagator: for V/j = R d 1 x [U,tf] and appropriate boundary 
conditions at spatial infinity, 

W[(<p f ,n),V fi ] = (<p f \e- H ( t f-W h \<p i ). 

(A3) W[(p, V] is translation and rotation invariant: i.e. under an isometry / : R d — > R d , 

W[<pof-\f(V)] = W[<p,V]. 

(A4) There is a functional derivative — on FCE) whose action on W[ip, V] can be approximated as follows: 

dip 

m 



^ d _, d n W a a ^o f , v , > S n W 



To evaluate the path integral (36), it is useful to arrange the field variables from each point in vectors <j) 
and write the action as 

S[cj>,V a } = ^<j>-B a( f> + c a -(f> + d a . (37) 

The boundary fields (p a are contained in the vectors c a and d a respectively. The action is bounded from 
below, so for each <p a , there is at least one solution cj) c i of the Euclidean equations of motion 

dS 

Qjlfa, V a ] = B a </> + c a = . 

If B a is non-degenerate, the solution is unique and one can define the Hamilton function 

S[<p a ,V a ] := S[4> cl ,V a ] 
for the discrete Euclidean system. We assume, in fact, that 
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(A5) The matrix B a is non-degenerate and the Hamilton function S[tp a , V a \ is analytic in ip a . 



The change of variables £ := <fi — <fi c i renders the action (37) quadratic: 

S[Z,V a ] = ±Z-B a Z + S[<p a ,V a ] 

The integral (36) becomes Gaussian and gives 



W a [<p a ,V a ] = 




Therefore, by (A5), the regularized kernel W a [p a , V a ] must be analytic in ip a , which will be used in section 
5.1 when deriving the Schrodinger equation. 

Remark: In (Al) and (A5) we have formulated the continuum limit in terms of pointwise convergence, 
i.e. by separate convergence for each boundary field <p in F(T,). According to (38), the field dependence 
of the regularized kernel resides only in the Hamilton function S[tp a , V a ]. The latter converges against the 
continuum function S[(p, V], which is defined pointwise. Thus, it is plausible to assume that the continuum 
propagator W[<p, V], too, is a pointwise function on F(E). When further developing the formalism, pointwise 
convergence is likely to be replaced by convergence in a Hilbert space norm or other measures which only 
distinguish between equivalence classes of boundary fields. For the purpose of this article, however, it is 
sufficient and simplifies notation if we use convergence on single fields. 



5 Generalized Schrodinger Equation 

The propagator W depends on a spacetime region V and a field (p specified on the boundary S. Thus, as 
for the Hamilton function in section 3, one can define a deformation derivative: using the same notation as 
there, we set 

L N W lv ,V]:=^ W ^- W ^ V l 

s^O S 

In this section, we derive that 

hL N W[<p,V\ = (-H N [ip,h^,V] + P N [<p,h^,V])w[<p,V], (39) 

where 

H N [<p,*,V] := ^cEA^-^r 2 + ^(Vs^) 2 + ^mV) , 
P N [<p,n,V] := -^cEiVy ■ Ve^tt. 

When V = R d_1 x [U,tf] and N±_\ tf = 1, N±\ ti = 0, ATy = 0, this yields the ordinary Euclidean Schrodinger 
equation 

The strategy of the derivation: using assumption (A5) and rotation invariance (A3), we show that the 
regularized propagator satisfies a lattice version of equation (39) when V is deformed along flat parts of its 
boundary. The central step is analogous to the calculation Feynman used when deriving the Schrodinger 
equation from the path integral of a point particle [14] (see also chap. 4, [17]). Due to (Al) and (A4), the 



16 



<> o 9 o- 



O — i — 6 — i — 6 — i — 6- 



n 



Figure 9: Addition of a single layer. 

discrete equations have the continuum limit (39). To cover also the case, when deformations are applied to 
curved sections of E, we approximate E by a triangulation, apply (39) to each triangle and let the fineteness 
of the triangulation go to zero. 

For simplicity, the argument is formulated for bounded volumes below. The generalization to infinitely 
extended V is straightforward. 



5.1 Discrete Schrodinger Equation 

Consider a lattice diagram in which part of E a coincides with a hypersurface H a of the lattice L a . Let 
n denote the normal vector of H a . The simplest way of modifying such a diagram is to add a {d — 1)- 
dimensional layer of points along H a n E a (see Fig. 9). The old boundary points adjacent to the layer 
become interior points. We describe this by a lapse function N a : E a — > {0, 1} which indicates for any given 
point of the boundary if a new point will be linked to it or not. Then, the function 



o-„ 



L a , x i— > x + aN a (x) n 



is the discrete flow associated to the deformation of the boundary E a . Define the new diagram and its 
boundary by 

V^.= V a Ua a (Z a ), E' a := a a (E a ) . 

The set 

Ex ^^(l) 

is the part of E a which is moved and becomes E' x := <r a (Ei), while 

So ^AT-^O) 



remains unchanged. As in the continuous case, we choose the new boundary field to be the pull-forward of 
the old one, that is, 

The resulting path integral is 



° 0- a 1 



w a y a x] = [ n v 7 ^) f n fcad^x) 



l d 

xexp<;- ¥ 



- ^s[<p, V a 
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(Recall that IQZ'i) is the set of links between points of E'i_.) The same can also be written as a convolution 
with the original propagator: 



J] I C a d4>(x))exp\~ 



^ 1 ( ip a (x) - <j){x) 



+ E lwi<p*) 2 + E \ m2 rti? 



W a [&a\i: A),Va] 



Following Feynman's derivation of the Schrodinger equation, we introduce new variables 



7 d-2 



ft. 



■(</>(x) - <A*(z)) , x G Si 



and get 



WaWaX] 



--1/V27T 



n / c -vQ^ exp (4 s ^ - t f s w + e ^v«(x)] } 



.xeEi 
x W a 



Next we apply Laplace's method to obtain an asymptotic expansion of this expression (see e.g. chap. 11, 
[17]): the dominant contribution to the Gaussian integral comes from an a-dependent interval [— e a ,e a ]' Sl ' 
around £ = 0. (|£i| denotes the number of points in Si.) The integral outside is exponentially damped 
for a — > and neglected. Within the interval, one can Taylor expand W a in £ and reverse the order 
of integration and Taylor expansion. To evaluate the integration for each term, the integration range is 
extended back to its full size: this introduces an error in each term of the sum and convergence is lost, but 
the expansion is still valid asymptotically for a — ► 0. 

For the integrations and estimates, we use the following formulas: 



J dyy n e-y 



-y 2 /2 _ 



(n — 1) (n — 3) • • • 3 • 1 ■ \^2tt , n > and even , 
, n odd , 



J dyy n e~ y2/2 = O^ 1 e" £ ' /2 ) as a -► . 



±£a 

We set e a = 1/a. Consider first the integral outside the chosen interval: 



ni / -p Ue 

^ V H\[-l/a,l/«] / V 



— 6XP " ft 



(40) 



(41) 



E ^( V ^a) 2 + E \™ 2 &(?) 

«e/(Ei) xeEi 



X W a 



a|s , <Aj|Ei + \/V a<i 2 Cj i 



(42) 



The second exponent vanishes in the continuum limit. For W a , we employ formula (38) and replace 
exp(— S[. . . , V a ]) by 1, as the action is positive. The determinant and C a -factors are together of order 
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0(1), since, by assumption, (38) approaches a finite continuum limit. Thus, the modulus of (42) is smaller 
than 

(41) _ (\ _|Sl| X 



0{l) < O -e 2? asa^O. 



^esx ■ RM _ 1/a;1/a] 
In the integral over [— I /a, l/a]l El l, we Taylor expand W a in £: 



1 



M^i ^ J- 1 /* 



d£(x) | exp | -- 2^ — | exp 



2 " a 



a 

7T 



E \(ywa?+ E H™Va(*) 



ft 



By assumption (A5), W a is analytic in the field variable, so the Taylor expansion converges uniformly and 
we are allowed to integrate each term of the series separately. We also set the limits of integration back to 
plus and minus infinity. This does not affect the asymptotic property of the series, since for each term the 
resulting error is only exponentially small: for example , the linear term yields 



1 



1/a 



dttx) e 



E 



d-l 



dW, 



9(a rf -Va(x)) V ad ~ 2 



O 



- i 

: e 2^ 



Va d ~ 2 



because of (41) and (A4). Then, we can use equation (40) to do the Gaussian integration in each term of 
the asymptotic series. Each integration, that is, each point x € Si, leaves an overall factor y/2-ir. Terms 
with an uneven number of £ variables (of the same point) vanish. We obtain 



w a y a x\ ~ ^p 



E ^(v^) 2 +E^V^)l} 



x I W a [ Va , V a ] + E 2 ^2 • "^2 + E E c(n) 



<9y? a (x) 2n (a' 



2n f7,<2-2^ra 



where the c(n)'s are numerical coefficients. If we write ip a as the pull-back ip' a o cr a = cr*^ of and use 
also the lapse function N a , the final result becomes 

w a y a x\ ~ exp(-i e « d_1 ^4(v^X) 2 + E « d_1 « iv 4 m2(<T ^« )2(x) 1 

[ L«e/(s a ) zes a z J J 



+ E 



d(n-l) 



c(n) 



n=2 



a(a rf -Va(^)) 



2 5(a^Va(x)) 2 



2n 



(43) 



Iteration 

Suppose now that the deformed set V a does not arise from the addition of a single layer, but from a contin- 
uous deformation V s of the orginial volume V. That is, we want to compare W [^„, V£] against W a [ip a , V a \. 
To make the calculation tractable, we require that the vector field N vanishes outside a neighbourhood U of 
a boundary point x £ S, and that within U the boundary £ is flat. Denote this part of £ by £[/ := S PI U. 
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Figure 10: Diagram for V£. 



By rotation and translation invariance ((A3)), we can orient V such that E[/ coincides with a hyperplane 
of the lattice L a . Let us begin by considering the case where the lapse N± is positive, that is, V C V s . For 
small enough s, the typical diagram for W [</?„, V£] looks like Fig. 10 (or its higher-dimensional equivalent), 
where along the normal direction n each point of E a is in one-to-one correspondence with a point of E*. 
(Note that in the limit s — ► 0, the slope of E s against E becomes arbitrarily small.) The new boundary 
E* can be built from E a by repeatedly adding single layers as described previously. Thus, we can iterate 
formula (43) to obtain a relation between IKbaiKl an d ^aka^ai ^o] where now, a a is the concatenation 
of all single-step flows. When collecting the various terms of the iteration, the lapse functions for each step 
add up to the total lapse function N a . We order the result in powers of aN a and a: 



Wa[<p s a ,V°] = W a [a* a <p s a , V a \ + Yl a d ~ l aN a 



n 



d 2 W a 



[<<P 8 a,V a ] 



1 

h 



2 d{ a d-^ a { x )Y 



lel(E a ) xeT, a 
+ 0(a d (aN a )) + 0({aN a ) 2 ) as a, s -► . 

Note that the displacement vector aN a approaches sN± when both a and s become small, i.e. 

aN a = sN ± + 0(s 2 ) + O(a) . 

If N is normal to Ej/, the discrete flow cr a approximates the continuous one, ct s , and 

= (f a + 0(a) . 
In general, iV has also a tangential component, so 

<T* a <p' a = <p a - sNtfV^a + 0{S 2 ) + O(o) , 



(44) 



(45) 



(46) 



as can be seen from the arrow diagram in Fig. 10. Plugging (45) and (46) into (44), one arrives at a 
regularized form of the Euclidean Schrodinger equation: 



W a [<p S a , V*] - W a [(f a 



,v a ] 



O a W a [cp a , V a ] + O(s) + 0{a) + 0{a/s) , 



(47) 
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Figure 11: Lapse with positive and negative sign. 



where O a is the operator 

~\ E a^i^s)^,^) 2 (48) 

An analogous argument applies to the case of negative lapse N±. For mixed diagrams as in Fig. 11, both 
types of calculations can be combined to give (47) for lapses of arbitrary sign. 



6. =- -i £ 



5.2 Continuous Schrodinger Equation 



Choose N as before, i.e. with support on a neighbourhood U of x G S where £ n ?7 is flat. We want to 
show that 

L N W[<p, V] = lim = 6w[ v] (49) 

s^O S 



where 



O 



N ± (x) 



2 8 V W + 2^ )2{x) + 2 m ^ {x) 



+ N\\ (x) ■ Ve¥>(s) ft 



5 



5v?(x) 



Stated more explicitly, this means that for any e > there is an sq > such that 

W[ip s ,V s ] - W[<p,V] 



OW[<p, V] 



e 



for all s < sq 



(50) 



To obtain an upper estimate on the left-hand side, we insert regularized propagators and operators in a 
suitable way, and then apply the triangle inequality: 



LHS of (50) = 



-(w[p a , V s ] - W a [<p' a , V a s ] + W a [<p' a , V°] - W a [<p a , V a ] + W a [<p a , V a ] - W[<p, V]) 



OaWa^a, V a ] + O a W a [<p a , V a ] - OW[ip, V] 



< 



S 
+ 

+ 



W[<p s ,V s ]-W a [<p s a ,V a s ] + 

dw[<p,v\-d a w a [<p a ,v a ] 

W a [<p s a ,V a s ]-W a [<p a ,V a 



W[<p,V]-W a [<p a ,V a } 



- d a w a [<p a ,v a ] 
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By assumption (Al) (existence of the continuum limit), the first two terms become smaller than e/4 when 
the lattice constant a is smaller than some a s > 0. The partial derivatives and potential terms in (48) 
approach their continuum analogues as a — > 0, so there is also an ao > such that 



< - for all a < ao ■ 



OW[<p,V]-O a W a [<p a ,V a ] 

The regularized Schrodinger equation tells us that for s smaller than some sq, there is an a' s > such that 
W a [^Vi\-W a [ip a ,V a ] 



- O a W a [ip a ,V a ] 



< 



for all a < a' . 



Thus, for any s < sq, we can choose a < min{a s , ao, a' s } and the left-hand side of (50) must be smaller than 
e. □ 



5.3 Curved Boundaries 

As it is based on the lattice equation (47), the previous derivation applies only when flat sections of the 
boundary £ are deformed. We do not know how to extend the lattice calculation to the case where both 
initial and deformed surface are curved. Below we give an argument which circumvents this difficulty, but 
requires additional assumptions. The idea is to approximate the curved boundary by a triangulation, apply 
the variation to each of the flat triangles and add up the contributions. 

Let Ts be a triangulation of £ with fineness <5: that is, when two 0-simplices are connected by a 1- 
simplex, their metric distance is at most S. Let {£<*} denote the set of (d — l)-simplices S a C S of the 
triangulation. The corner points of each such simplex S a define a (d — l)-simplex in R d which we call 
£a„- The hypersurface £a := U q £a q approximates S and encloses the volume Va- We can view V as a 
deformation of Va and find a flow 

^RxK^K, (t, x) i-> p(t, x) = p t (x) 

such that Pi(Va) = V and Pi(Sa q ) = S a . We equip Sa with the boundary field <p& := p\ip = ip o p l5 the 
pull-back of (p under this flow. Motivated by equation (49) for flat surfaces, we assume that the difference 
between W[ip, V] and VF[<£>a, Va] is of the order of the volume difference between V and Va- 

W[ip,V] = W[puVA,Pi(V A )} 

= W[<pa,V a ]+0(\V-Va\). (51) 

Next we introduce "characteristic" functions Xa '■ ^ d —> K- with the property that 

Xa(x) = 1 for x G Sa q , 
Xa(x) =0 for x G S A(3 , a / (3 , 
and J2a Xa(x) = 1 for all x G R d . 

Using these functions, we can decompose the deformation field ./V according to 

a a 

Each component N a is a discontinuous vector field and gives rise to a discontinuuous flow within 1R, d . 
Suppose that by a suitable limiting procedure, one can define L^ a such that equation (49) holds and 

a 
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Then, equation (51) becomes 



L N W[p, V] = J2 L Na W[ip A , V A ] + 0(\V -V A \). 



By construction, the vector fields N a are only nonzero on the flat simplices £ a . Therefore, our result for 
flat surfaces (equation (49)) is applicable and yields 



L N W[<p, V] 



+ 0(\V-V A \) 



K z 5 2 1 



2 5ip 



2 + ^(Vs^a) 2 + 2 m Vi )+ N W V £^A fliz 



5_ 

5ip 



W&a,Va] 



h 2 5 2 



N ± - 



1 



S 



2 V + 2 (V ^ )2 + 2 mV +Ar ll" V ^% 



w>,y] + o(|y-u A |) 



In the (5^0 limit, |V — Va| goes to zero and one recovers the generalized Schrodinger equation for curved 
boundaries. 



6 Summary 

We have proposed an exact definition for a Euclidean free scalar propagator W[ip, V] which "evolves" 
wavefunctionals of fields along general spacetime domains V . Our main result is a derivation of the evolution 
equation 

hL N w[<p,v\ = [-H N [^n^v] + p N [ v ,hj-,v]^w[p,v}. (52) 

This equation describes how W[ip, V] varies under infinitesimal deformations of V generated by a vector 
field N. The variation is given by the action of two operators: one is related to the field Hamiltonian and 
arises from normal deformations of the boundary S = dV. The second operator results from tangential 
deformations and generalizes the field momentum. 

We showed also that the Hamilton function of the classical system satisfies an analogous Hamilton- Jacobi 
equation 

L N S[<p, V] = H N [<p, ^,V]+P N [ V ,^,V}. (53) 

When the boundary S consists of two infinite hyperplanes at fixed times, (52) and (53) reduce to the 
standard Schrodinger and Hamilton-Jacobi equation in their Euclidean form. 

The derivation of eq. (52) is based on assumptions which we consider plausible, but are not proven. 
Most importantly, we have not shown that the proposed regularization of the propagator has a well-defined 
continuum limit. A description for converting the Euclidean to a Lorentzian propagator is missing. As 
described in section 2, we expect that an evolution equation analogous to (52) holds also for Lorentzian 
propagators. We emphasize that such state evolution may, in general, be non-unitary and nevertheless 
admit a physical interpretation. 
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A Local Form of Hamilton-Jacobi and Schrodinger Equation 



In the text we have presented the generalized Hamilton-Jacobi and Schrodinger equation in an integral 
form. They can be also expressed locally, and below we explain how the two representations are related. 
The local notation is used in [3] and [7]. 

Both the Hamilton function S and the propagator W depend on the volume V. The latter is enclosed 
by the boundary S. Consider a parametrization of S, i.e. a map 

x : P — > S , t i ^ x(t) 

from a (d — l)-dimensional manifold P to S. Provided it is clear on which "side" of S the volume lies, one 
can view S and W as functions of S, or equivalently, as functionals of the parametrizing map r x(r). 
The other variable of S and W is the boundary field <p : £ — > R: we may replace it by its pull-back ^ to P, 
so that S 1 and W are completely expressed in terms of quantities on the parameter manifold P: 

ip(r) = <p(x(t)) , t e P, 

S = S[<p(t),x(t)], W = W[<p(t),x{t)]. 

In section 3, we defined the deformation derivative Ljv which acts by infinitesimal diffeomorphisms and 
pull- forwards of V and ip respectively. A moment's reflection shows that in the new notation the same 
effect is achieved by applying a variation 

Sx(t) = sN(x(t)) 
to the function x(t) while leaving ip(r) fixed. Therefore, 

Our explicit result for the Hamilton-Jacobi equation (p. 10, eq. (22)) can be rewritten as 



(54) 



L N S[<p, V] 



d^TN^xir^ln^xir)) 



V M ^(r) 



1 / SS 



1 



1 



2VMr); +2^) 2 + 2^ 2(T) + ^ (r)) 
5S ) 



(55) 



where on the right-hand side S is a functional of the new variables. Comparison with (54) gives the equation 



5S 



8x»(t) 



n ^{x{r)) 



1 



5S 



1 



1 



5S 



(56) 



It describes how S behaves under local variations of the boundary S. By the same reasoning, we arrive at 
a local Schrodinger equation for the kernel W: 



5W 
5x^{t) 



h 2 5 2 W l.„ ± . 1 2j2/ ' 

+ 2 (V * (T)) + 2™ * (T) 



5W 

5<p(t) 



(57) 
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